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ABSTRACT 

We reanalyze the cosmological constraints on the existence of a net universal 
lepton asymmetry and neutrino degeneracy. We show that neutrinos can begin 
to decouple at higher temperatures than previous estimates due to several 
corrections which diminish the neutrino reaction rate. For sufficiently large 
degeneracy, neutrino decoupling can occur before various particles annihilate 
and even before the QCD phase transition. These decoupled neutrinos are 
therefore not heated as the particle degrees of freedom change. The resultant 
ratio of the relic neutrino-to-photon temperatures after annihilation can 
then be significantly reduced by more than a factor of two from that of the 
standard nondegenerate ratio. This changes the expansion rate and subsequent 
primordial nucleosynthesis, photon decoupling, and structure formation. In 
particular we analyze physically plausible lepton-asymmetric models with large 

and v T degeneracies together with a moderate v e degeneracy. We show that 
the nucleosynthesis by itself permits very large neutrino degeneracies < £j, , 
iv T < 40, < £„ B < 1.4 together with large baryon densities 0.1 < ^lbh\ Q < 1 
as long as some destruction of primordial lithium has occurred. We also show 
that structure formation and the power spectrum of the cosmic microwave 
background allows for the possibility of an Q = 1, Q\ = 0.4, cosmological model 
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for which there is both significant lepton asymmetry (\£ v \ = \£ Vr \ ~ 11) and a 
relatively large baryon density (O&/t§ ~ 0.2). Our best-fit neutrino-degenerate, 
high-baryon-content models are mainly distinguished by a suppression of the 
second peak in the microwave background power spectrum. This is consistent 
with the recent high resolution data from BOOMERANG and MAXIMA-1. 



Subject headings: cosmology: lepton asymmetry, J7 b - cosmic microwave 
background 



1. INTRODUCTION 



The universe appears to be charge neutral to very high precision QLyttleton fc Bond! 



|1959|) . Hence, any universal net lepton number beyond the net baryon number must reside 
entirely within the neutrino sector. Since the present relic neutrino number asymmetry is 
not directly observable there is no firm experimental basis for postulating that the lepton 
number for each species is zero. It is therefore possible for the individual lepton numbers L\ 
of the universe to be large compared to the baryon number of the universe, B, while the net 
total lepton number is small (L ~ B). Furthermore, it has been suggested that even if the 
baryon number asymmetry is small, the total lepton number could be large in the context 
of SU(5) and SO (10) Grand Unified Theories (GUT's) ( |Harvey fc Kolb l98l| ; [Fry fc Hogari 



1982| ; polgov fc Kirilova 1991| ; polgov 1992| ). It has also been proposed recently flCasas. 



Cheng, fc Gelmini 1999|) that models based upon the Affleck-Dine scenario of baryogenesis 



( Affleck fc Dine 1985Q might generate naturally lepton number asymmetry which is seven 



to ten orders of magnitude larger than the baryon number asymmetry. Neutrinos with 
large lepton asymmetry and masses ~ 0.07 eV might even explain the existence of cosmic 
rays with energies in excess of the Greisen-Zatsepin-Kuzmin cutoff QGrciscn 1966| ; Gelmini 



fc Kusenko 1999|) . Degenerate, massive (2.4 eV) neutrinos may even be required ( [Larscn 



fc Madsen 1995| ) to provide a good fit to the power spectrum of large scale structure in 
mixed dark matter models. It is, therefore, equally important for both particle physics and 
cosmology to carefully scrutinize the limits which cosmology places on the allowed range of 
both the lepton and baryon asymmetries. 

The consequences of a large lepton asymmetry and associated neutrino degeneracy 
for big bang nucleosynthesis (BBN) have been considered in many papers. Models with 
degenerate v e ( |Wagoner, Fowler, fc Hoyle 1967| ; [lerasawa fc Sato 1985| ; |Scherrer 1983| ; |1988 



Kajino fc Orito 1998a| ), both v e and ( |Yahil fc Beaudet 1976| ; |Bcaudet fc Goret 197£ 



Bcaudet fc Yahil 1977| ), and for three degenerate neutrino species ( pavid fc Reeves 1980| ; 
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Ulive et al. 19911 ; |Kang fc Steigman 199^ ; [Starkman 1992] ; |Ka.jino fc Orito 1998b| ; |Kim fc 
Lee 1995| [Kim, Kim fc Lee 1998| ) have been analyzed. The effects of the degeneracy of 
electron type neutrinos on inhomogeneous BBN models were also considered in Kajino & 
Orito (1998a, 1998b). Constraints on lepton asymmetry also arise from the requirement 
that sufficient structure develop by the present time QKang fc Steigman 1992|) and from 
the power spectrum of fluctuations in galaxies ( |Larsen fc Madsen T995|) and the cosmic 
microwave background temperature (|Kinney fc Riotto 1999| ; [Lesgourgues fc Pastor 199"9 



Hannestad 20001) . 



The present work differs from those listed above primarily in that we carefully 
reexamine the neutrino decoupling temperature. Although neutrino degeneracy has been 
treated in considerable detail previously (cf. [Kang fc Steigman 1992| ) we reexamine several 
issues which were not treated in that paper or elsewhere. We show that there are corrections 
to the neutrino reaction rates which cause the neutrino decoupling temperature to be 
significantly higher than previous estimates. These corrections include a proper accounting 
of effects of degeneracy on the abundances of the weakly- reacting species present. We 
also consider corrections in the the electron and photon effective masses and the relative 
velocities of interacting species. In addition we make a more careful continuous treatment of 
the various particle annihilation epochs. These corrections are important since it becomes 
easier for neutrinos to decouple before the annihilation of various particles and even before 
the QCD transition. The fact that neutrinos may decouple when there were many particle 
degrees of freedom causes the relic neutrino temperature to be much lower than that of 
the standard nondegenerate big bang. This allows for new solutions for BBN in which 
significant baryon density is possible without violating the various abundance constraints. 
We find that for decoupling temperatures just above the QCD epoch it is also possible to 
find models in which the structure constraint and even the CMB power spectrum constraint 
are marginally satisfied. 

In this paper we explore a physically plausible range of baryon and lepton asymmetric 
cosmological models in which the electron neutrino is only slightly degenerate while 
the muon and tau neutrinos have nearly equal (but opposite-sign) significant chemical 
potentials. We use these models to deduce new cosmological constraints on the baryon 
and lepton content of the universe. We emphasize that previous studies of BBN and the 
CMB temperature fluctuations with highly degenerate neutrinos have not exhaustively 
scrutinized the possibility of important effects from the lower implied neutrino temperature 
when neutrinos decouple before various particle annihilations and/or the QCD epoch. 



In what follows we first discuss here how the observable BBN yields in a neutrino- 
degenerate universe impose bounds on the baryon and lepton asymmetries which still 
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allow a large neutrino degeneracy and baryon density (even f2 h = 1). We show that there 
exists a one parameter family of allowed neutrino-degenerate models which satisfy all 
of the constraints from BBN. We next examine other cosmological non-nucleosynthesis 
constraints, i.e. the cosmic microwave background (CMB) fluctuations and the time scale 
for the development of structure. We show that these constraints can be marginally satisfied 
for a limited range of highly degenerate models from the fact the relic neutrino temperature 
is much lower than that of the standard nondegenerate big bang. We also discuss how a 
determination of the neutrino degeneracy parameters could constrain the neutrino mass 
spectrum from the implied neutrino contribution Q u to the closure density. 



2. LEPTON ASYMMETRY & NEUTRINO DECOUPLING 

In this section, we review the basic relations which define the magnitude of neutrino 
degeneracy and summarize the cosmological implications. Radiation and relativistic 
particles dominate the evolution of the early hot big bang. In particular, relativistic 
neutrinos with masses less than the neutrino decoupling temperature, m v < 0(Tb) ~ MeV, 
contributed an energy density greater than that due to photons and charged leptons. 
Therefore, a small modification of neutrino properties can significantly change the expansion 
rate of the universe. The energy density of massive degenerate neutrinos and antineutrinos 
for each species is described by the usual Fermi-Dirac distribution functions f u and f D , 

Pv + Pp = 

1 r°° i 

— j o dp p 2 ^p 2 + ml(U(p) + U(p)). (1) 

where, p denotes the magnitude of the 3-momentum, and m v is the neutrino mass. Here and 
throughout the paper we use natural units( H — c — Jzb — 1). The distribution functions are 

Hp) = ' 



exp(p/T v -£ v ) + l ^ 



Hp) = ' 



exp (p/T u + £,) + 1 ' 



where the degeneracy parameter £„ is defined in term of the neutrino chemical potential, 
fi u , as £j, = p v jT v . It will have a nonzero value if a lepton asymmetry exists. In the limit of 
massless neutrinos, the energy density in neutrinos becomes 



2 



7 7T , - 



15 (i v .\ 4 30 /f„. 

1 + y(v) +t(v 



(3) 



- 5 - 



from which it is clear that a neutrino degeneracy in any species tends to increase the energy 
density. 

The net lepton asymmetry L of the universe can be expressed to high accuracy as 

l=e,n,T (A\ 

n Vl — riv l 
L>\ — • 

This is analogous to the baryon-to-photon ratio 77 = (tie — ri§)/n 7 . Here, n v% (n p J are 
the number densities for each neutrino (anti-neutrino) species, n 7 is the photon number 
density, and % (n B ) is the (anti) baryon number density. Once the temperature drops 
sufficiently below the muon rest mass, say T < 10 MeV, all charged leptons except for 
electrons and positrons will have decayed away. Overall charge neutrality then demands 
that the difference between the number densities of electrons and positrons equal the proton 
number density. Hence, any electron degeneracy is negligibly small and any significant 
lepton asymmetry must reside in the neutrino sector. After the epoch of annihilation, 
the magnitudes of the lepton and baryon asymmetries are conserved. They are equal to 
the present value in the absence of any subsequent baryon and/or lepton number- violating 
processes. 

Elastic scattering reactions, such as i>\ {pi) + P 1 v% {pi) + 1"^, keep the neutrinos in 
kinetic equilibrium. Annihilation and creation processes which can change their number 
density, like v\ + f?i <-> I + I, vi + V <r-> vy + I, etc, maintain the neutrinos in chemical 
equilibrium. When the rates for these weak interactions become slower than the Hubble 
expansion rate, neutrinos decouple and begin a "free expansion". Their number densities 
continue to diminish as 1/R 3 and their momenta red-shift by a factor 1/R, where R 
is the cosmic scale factor. However, this decoupling has no effect on the shape of the 
distribution functions. Relativistic neutrinos and antineutrinos continue to be described by 
the Fermi-Dirac distributions of Eq. (fj). Since the individual lepton number is believed to 
be conserved, the degeneracy parameters £^ remain constant after decoupling. 



2.1. Neutrino Decoupling Temperature 

When one estimates the present density of relic neutrinos one must consider the effect 
of the changing number of degrees of freedom for the remaining interacting particles. For 
example, once the neutrinos are totally decoupled, they are not heated by subsequent pair 
annihilations. Hence, their temperature T v is lower than the temperature T 7 of photons (and 
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any other electromagnetically interacting particles) by a factor y u = T u /Tj. In the standard 
non-degenerate cosmology, with three flavors of massless, non-degenerate neutrinos which 
decouple just before the e ± pair annihilation epoch, the present ratio of the neutrino to 
photon temperatures is given by y u = (4/1 l) 1 / 3 . 

Neutrinos and anti-neutrinos drop out of thermal equilibrium with the background 
thermal plasma when the weak reaction rate per particle, T, becomes slower than the 
expansion rate, H. Since the decoupling occurs quickly, it is a good approximation to 
estimate a characteristic decoupling temperature To, at which the slowest weak-equilibrating 
reaction per particle T becomes slower than the expansion rate. 

The expansion rate is simply given by the Friedmann equation 



H=M/3)*G'£p i , (5) 

where the sum is over all species contributing mass energy (mostly relativistic particles). 
Obviously, the universal expansion in neutrino-degenerate models is more rapid because of 
the higher neutrino mass-energy density (e.g. Equation |3[). This causes the decoupling to 
occur sooner and at a higher temperature than in the non-degenerate case. 

Our estimates of the weak decoupling temperature as a function of neutrino degeneracy 
are shown in Figure 1 and compared with the previous estimates Kang & Steigman (1992). 
Our decoupling temperatures are significantly higher than those estimated in Kang & 
Steigman (1992) for reasons which we now explain. 

Significant neutrino degeneracy will cause the weak reaction rate per particle, T, to be 
slower than the nondegenerate case because the neutrino final states are occupied ([Kang] 



|fc Steigman 1992] ; |Freese, Kolb, fc Turner 19831 ). More importantly, the existence of a 



lepton asymmetry represented by a finite positive neutrino chemical potential will cause the 
number density of antineutrinos to be significantly less than that of neutrinos. Hence, the 
annihilation rate per neutrino T v = np(a u p) can be significantly less than the reaction rate 
per antineutrino, electron or positron. 

In Kang & Steigman (1992) the weak reaction rate was estimated from the reaction 
rate per electron annihilation, e + + e~ — > vi + z7j. We will call this T e . One could however 
also consider the reaction rate per positron, T e +, the reaction rate per neutrino species i, 
T u ., or the reaction rate per antineutrino, T p .. In the limit of high temperature and no 
neutrino degeneracy, all of these reaction rates are identical. However, when for example 
the neutrinos have a large positive chemical potential, the neutrino annihilation rate 

— n u i (o'v) becomes the rate determining reaction. This is because of the scarcity of 
anti- neutrinos with which to interact. It is then this reaction rate and not r„ which should 
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be used to determine when the weak reactions are no longer able to maintain chemical 
equilibrium. 

In Kang k, Steigman (1992) expressions were derived for the phase space integrals for 
(av) which reduce to 

where is a phase space integration factor given in that paper, ( is the Riemann zeta 
function, and Gp is the Fermi coupling constant. The quantities a and b relate to the 
Weinberg angle, a = 1/2 + 2sin 2 #^ for electron neutrinos, and a = 1/2 — 2sin 2 #v^ for 
muon and tau neutrinos, while 6 = 1/2. 

Using this, we derive the following expression for the decoupling temperature 
corresponding to each neutrino species, i, 



4.38 lV 3 r V^//(0 



a 2 + b 2 



1/3 / 

x U) MeV , (7, 



where % denotes the matter (or antimatter) counter part for the species under consideration. 
The formula A. 8 of Kang & Steigman (1992) omits the latter density ratio factor which 
is only unity in the limit of no degeneracy. Thus, their T used to deduce the decoupling 
temperature was too large. We have checked the above result against a full eight-dimensional 
integration of the Boltzmann equation and find our results to be correct. 

There are also two other differences between our results and those of Kang & Steigman 
(1992): The first is that we have included finite temperature corrections to the mass of the 
electron and photon ([Fornengo, Kim, fc Song 1997|) . The second is that we have calculated 



the average neutrino annihilation rate in the cosmic comoving frame. In this frame the 
M0ller velocity is used instead of the relative velocity for the integration of the collision 
term in the Boltzmann equations flGondolo, fc Gelmini 199"T| ; |Knqvist et al. 1992]) . These 



corrections also slightly increase the neutrino decoupling temperatures. Therefore even in 
the nondegenerate case we find slightly higher decoupling temperatures than those of Kang 
& Steigman (1992): 

T D (f„ = 0) ~ 2.33 MeV for u e , 

T D (£ V = 0) ~ 3.90 MeV for v^ T . (8) 

These values are in good agreement with those of Enqvist et al. (1992). If we remove the 
above corrections, we recover the decoupling temperatures of (Kang & Steigman 1992), i.e. 
T D {£ V = 0) ~ 1.98 MeV for u e , and T D (£ U = 0) ~ 3.30 MeV for v^ T . 
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2.2. Relic Neutrino Temperature 

One does not need to increase the decoupling temperature by much before photon 
heating by annihilations becomes important for determining the relic neutrino temperature. 
For illustration, Figure shows the ratio of muon to photon energy densities as a function 
of temperature in units of the muon rest mass m M = 105 MeV. A similar curve could be 
drawn for any massive species. One can see that even at a temperature of only 20% of the 
muon rest mass, muons still contribute about 10% of the mass energy density, and hence, 
can affect the ratio of the photon to neutrino temperatures as these remaining muons 
annihilate. Combining Figures [I] and |2|, one can see that even for a degeneracy parameter 
of £„ ~ 6, the decoupling temperature is at 20% of m^. For the case of highly degenerate 
neutrinos > 10.8 and £„ T > 9.8), Td(^J can exceed the muon rest energy and even the 
QCD phase transition temperature. 

If the neutrinos decouple early, they are not heated as the number of particle degrees 
of freedom change. Hence, the ratio of the neutrino to photon temperatures, T u /T^, is 
reduced. The computation of the ratio of the final present neutrino temperature to the 
photon temperature is straightforward. Basically, since the universe is a closed system, the 
relativistic entropy is conserved within a comoving volume. That is; 

R 3 s = Constant , (9) 

where the entropy density s is defined, 

S = 2s f - ~^9eff-L , (10) 

and the sum is over all species present. Since the mass-energy is dominated by relativistic 
particles s can be written in terms of the effective number of particle degrees of freedom, 

9eff= E 9i(j) +1 E 9i(j) • (11) 

i=bosons i=fermions 

As each species annihilates, sR 3 remains constant. Therefore, the temperature of the 
remaining species increases by a factor of (^//^V^e// 61 ") 1 ^ 3 - This accounts for the usual 
heating of photons relative to neutrinos due to pair annihilations by a factor of (11/4) 1 / 3 . 
Note, that in the computation of g e ff from equations ( |T0"D and ( |TT]) it is important to 
evaluate the energy densities continuously (cf. Figure |^) and not simply assume abrupt 
annihilation as the temperature approaches the rest energy of each particle as was done in 
Kang & Steigman (1992). 

Figure [3] illustrates g e ff as a function of temperature from 1 MeV to 1 TeV. In each 
calculation g e ff depends upon the relic neutrino temperature and therefore when the 
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neutrinos decoupled. For purposes of illustration only, we drawn this figure for three 
neutrino types which do not decouple until ~ MeV. To construct this figure we have 
included: mesons (tt, p, <fi, u, r},rf)\ leptons (e, /i, r, z/ e , z/ T ); a QCD phase transition at 
150 MeV; light quarks and gluons (u,d, gluon); s-quarks (with m s = 120 MeV); c-quarks 
(with m c = 1200 MeV); b-quarks (with m b = 4250 MeV); W, Z-bosons (with m z = 80 
GeV); and t-quarks (with m t = 173 GeV). 

One can see from this figure that if neutrinos decouple before the QCD phase transition, 
there is substantial heating of photons due to the large change in degrees of freedom at this 
epoch. 

Figure § shows the final ratio of neutrino temperature today to that of the standard 
non-degenerate big bang for three neutrino flavors. The biggest drop in temperature for all 
three neutrino flavors occurs for £„ « 10. This corresponds to a decoupling temperature 
above the cosmic QCD phase transition. The low temperature is the result of the decrease 
in particle degrees of freedom during this phase transition. This discontinuity will have 
important consequences in the subsequent discussions. 

For comparison, the same computation from Figure 1 of Kang & Steigman (1992) is 
shown. One can see that there is a substantial difference between the two results. These 
differences stem from two effects. On the one hand the decoupling temperature is higher 
in the present work for a given £„. because of the various corrections as explained above. 
This is the reason that the discontinuity in relic neutrino temperature occurs at £ v T w 11 
instead of 15 as in Kang & Steigman (1992). The second difference is that here we have 
calculated the proper thermodynamic energy density for each fermion and boson species 
continuously (cf. Eq. 1), whereas in Kang & Steigman (1992) an approximation was made 
that each particle species instantaneously annihilates as the temperature drops below its 
rest energy. This is the reason for the step function appearance of the curve in Kang & 
Steigman (1992). 

Some discussion of book keeping here is in order as it explains why we observe the jump 
in relic temperature when neutrinos decouple above the QCD transition. Let us consider a 
neutrino which decouples just above the QCD transition. Then the relevant fermion degrees 
of freedom consist mainly of 3 q — q colors (treating the s-quark as relativistic) x 3 flavors 
contributing gf = 31.5, plus e + — e~ and fi + — /x~ plus 2 remaining nondecoupled flavors 
of v — v lepton pairs. These contributing another gf = 10.5. For the bosons there are 8 
gluons and 1 photon contributing a total g^ = 18. While the species are in equilibrium 
the total g e ff = 60. If no more species decouple before the electron pair annihilation 
epoch, the fermion degrees of freedom before electron-positron annihilation (at T ~ 1 
MeV) include only 2 flavors of v — v plus e + — e~ so that gf = 7, while for bosons only 



-10- 



photons remain so that gb — 2 and g e ff = 9. The relic neutrino temperature is therefore 
reduced by an additional factor due to these differences in particle degrees of freedom 
(T„/T 7 )/ [4/1 1] 1 / 3 = (9/60) 1 / 3 « 0.53 (cf. Figure §. 

For all three neutrino flavors the temperature begins to decrease relative to the 
standard value for a degeneracy parameter as small as £ ~ 5. This is because some relic 
\i — Jx pairs are still present even at temperatures well below the muon rest energy (cf. 
Figure |2|). The first neutrino species to be affected as the degeneracy increases are the 
and v T . They decouple at a higher temperature than v e even in the standard nondegenerate 
big bang because the electrons continue to experience charged- current interactions to 
lower temperature. The muon neutrinos exhibit a slightly different behavior than v T for 
degeneracy parameters £„ > 5 because the \i — Jjl density is large enough at the decoupling 
temperature for charged- current interactions to be significant. This maintains equilibrium 
for to somewhat lower temperatures even for degenerate neutrinos. This causes the 
decoupling temperature to be lower (cf. Figure P and the relic temperature to be slightly 
higher than the v T temperature for degeneracy parameters between 5 and 9. 



PRIMORDIAL NUCLEOSYNTHESIS 



3.1. Current Status 



Although the homogeneous BBN model has provided strong support for the standard 
hot big-bang cosmology, possible conflicts have emerged between the predictions of BBN 
abundances as the astronomical data have become more precise in recent years. One 
difficulty has been imposed by recent detections of a low deuterium abundance ( [Buries 
|fc Tytler M| |Burles k Tytler 1998 b| , see also [Levshakov, Tytler, k BurlesTDT^ ) in 
Lyman-a absorption systems along the line of sight to high red-shift quasars. The low D/H 
favors a high baryon content universe and a high primordial 4 He abundance, Y p > 0.244. 
This is inconsistent with at least some of the reported constraints from measurements of 
a low primordial abundance of 4 He, Y„ rs 0.235 ± 0.003, in low-metallicity extragalactic 



H II regions ( piive k Stcigman 1995 ; Stcigman 1996 ; Hata et al. 1996 ; piive, Stcigman 



fc Skillman T997| ; |Kajino k Orito 1998a| ; [Piembert k Piembert 2000]) . This situation is 
exacerbated by recent detailed analyses ( [Esposito et al. 1999| ; Lopez k Turner 1999) of the 
theoretical uncertainties in the weak interactions affecting the neutron to proton ratio at 
the onset of primordial nucleosynthesis. These results require a positive net correction to 
the theoretically determined 4 He mass fraction Y p of +0.004 to +0.005 or ~2%. We also 
note that the low deuterium abundance is marginally inconsistent with the 7 Li abundance 
inferred by measurements of lithium in Population II halo stars ( [Ryan et al. 19991 ; [Kajino 



- 11 - 



et al. 2000). Significant depletion of lithium from these stars, or a lower reaction rate for 
primordial lithium production may be required. 

Another potential difficulty has been imposed by recent X-ray observations of rich 
clusters QWhite et al. 19931 ; IWhite fc Fabian 1995fc Pavid, Jones, fc Forman 1995| ; |Bahcall, 



Lubin, fc Dorman 1995|). The implied baryonic contribution to the closure density is 



0.08 < f2 h h% 2 /Q M < 0.22 ( [lytler et al. 2000|) , where tt M is the total matter (dark plus 
visible) contribution, and /2.50 is the Hubble constant Ho in units of 50 km s _1 Mpc -1 . 
Consistency with the limits (0.03 < ^b^Jo < 0.06) from homogeneous BBN ( |Walker et 



|al. 1991| ; |Smith et al. 1993| |Copi, Schramm fc Turner 1995| |Schramm fc Mathews 1995 



TFT - 

Olive, Steigman, fc Walker 1999| ) then requires that 0.14 < Q^fh 50 < 0.75. Hence, matter 



dominated cosmological models (for example with Hq = 75 km s 1 Mpc 1 and Qm > 0.61) 
can be in conflict with BBN. 

Another possible conflict has emerged in recent measurements of the power spectrum 
of the cosmic microwave background (CMB) temperature fluctuations. A suppression of 
the second acoustic peak in the power spectrum has been recently reported in the balloon 
based CMB sky maps from the BOOMERANG ( [Bernardls, P. et al. 2000]) and MAXIMA- 1 
( [Hanany et al. 200U| ) collaborations. Indeed, the derivation of cosmological parameters 
from these new data sets ([Lange et al. 20001 ; |Balbi et al. 2000]) indicate a preference for a 



large baryon density. For example, the optimum cosmological models consistent with the 
BOOMERANG data (|Lange et al. 200q) indicate Q b h 2 = 0.032 ± 0.004 (Id) for fits in 
which there is no prior restriction on the range of Vt b h 2 . Similarly, the likelihood analysis 
based upon the MAXIMA-1 data flBalbi et al. 2000|) indicates tt h h 2 = 0.025 ± 0.005. These 



results are to be compared with the best current (la) limit from standard primordial 
nucleosynthesis without neutrino degeneracy, Q b h 2 = 0.019 ± 0.0024 ( |Qlive, Steigman,~fe] 



|Walker 1999| ; [l'ytler et al. 200U| ) which come from the deuterium abundance in Lyman-a 
clouds observed along the line of sight to background quasars. The CMB data indicate a 
baryon content which is at least 1-2 cr above the optimum BBN value. Indeed, Q b h 2 > 0.20 
would require a primordial helium abundance of Y p > 0.25 which is beyond even the most 
generous adopted limits from observation (|Qlive, Steigman, fc Walker 1999|) . Thus, it seems 



that that the CMB data hint at (though not necessarily prove) the need for a modification 
of BBN which allows higher values of Qbh 2 while still satisfying the constraints from light 
element abundances. Such conditions are easily satisfied by neutrino-degenerate models. 
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3.2. Neutrino-Degenerate BBN 



Most previous works have only considered the effects of neutrino degeneracy on the 
light elements 4 He, D, and 7 Li. Recently, the predicted abundance of other elements such 



as 6 Li, 9 Be, and n B in a neutrino-degenerate universe were also studied ( [Kim fc Lee 



|1995| ; [Kim, Kim fc Lee 1998| ). Here we investigate the effects of lepton asymmetry on the 
predicted abundances of heavier elements (12 < A < 18) as well as these light elements. 



Non-zero lepton numbers primarily affect nucleosynthesis in two ways ( |Wagone7 



bowler, fc Hoyle 1967fc [lerasawa fc Sato 1985fc [19881 [Scherrer 1981 ; |Yahil fc Beaudet 1976 
Beaudet fc Goret 197§ [Beaudet fc Yahil 1977| ; pavid fc Reeves 1980| ; piive et al. 1991 



Kang fc Steigman 1992| ; frjtarkman 199^ ; [Kajino fc Orito 19984 [Kim fc Lee Egg: [Rim 



Kim fc Lee 1998| ). First, degeneracy in any neutrino species leads to an increased universal 
expansion rate independently of the sign of (cf. Eqs. |] and |). As a result, the weak 
interactions that maintain neutrons and protons in statistical equilibrium decouple earlier. 
This effect alone would lead to an enhanced neutron-to-proton ratio at the onset of the 
nucleosynthesis epoch and increased 4 He production. 

Secondly, a non-zero electron neutrino degeneracy can directly affect the equilibrium 
n/p ratio at weak- reaction freeze out. The equilibrium n/p ratio is related to the electron 
neutrino degeneracy by n/p = exp{ — (AM/ T n ^ p ) — £„ e }, where AM is the neutron-proton 
mass difference and T n ^ p is the freeze-out temperature for the weak reactions converting 
protons to neutrons and vice versa. This effect leads to either increased or decreased 4 He 
production, depending upon the sign of L e or £ Ue . 

There is also a third effect which we emphasize in this paper. As discussed in the 
previous section, T v /T y can be reduced if the neutrinos decouple at an earlier epoch. This 
lower temperature reduces the energy density of the highly degenerate neutrinos during the 
BBN era, and hence, slows down the expansion of the universe. This leads to decreased 
4 He production. We show in the next section that the allowed values for ^ e ,^ , £ Ut and 
17b which satisfy the light-element abundance constraints are significantly changed for large 
degeneracy (£„ , £„ T > 10) compared to the results of previous studies. 



3.3. Summary of Light-Element Constraints 



The primordial light element abundances deduced from observations have been 
reviewed in a number of recent papers (|01ive, Steigman, fc Walker 1999[ |Nolett fc Burle; 



|2000| ; (Steigman 2000| ; [Tytler et al. 2000|) . There are several outstanding uncertainties. For 
primordial helium there is an uncertainty due to the fact that deduced abundances tend to 
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fall in one of two possible values, one high and the other low. Hence, for 4 He we adopt a 
wide range: 

0.226 < Y p < 0.247. 

For deuterium there is a similar possibility for either a high or low value. Here, however, 
we adopt the generally accepted low values of Tytler et al. (2000), 

2.9 x 10~ 5 < D/H < 4.0 x 10" 5 . 

For primordial lithium there is some uncertainty from the possibility that old halo stars 
may have gradually depleted their primordial lithium. Because of this possibility we adopt 
the somewhat conservative constraint: 

1.26 x 10~ 10 < 7 Li/H < 3.5 x 10~ 10 



3.4. Nucleosynthesis Model 

As we shall see, shifts in the relic neutrino temperature can dramatically affect the 
abundance yields ( [Kajino fc Orito 1998b|) . We now explore the parameter space of neutrino 



degeneracy and baryonic content to reinvestigate the range of models compatible with the 
constraints from light element abundances. 

For the present work we have applied a standard big bang code with all reactions 
updated up to A=18. [However, in the present discussion only reactions involving nuclei up 
to A=15 are significant.] In this way possible effects of lepton asymmetry on heavier element 
abundances could be analyzed along with the light elements. In this context a recent 
compilation of the nuclear reaction rates relevant to the production of n B ( Prito, Kajino 



|fc Oberhummcr 199S ) was useful because several important LiBeB(a,x) and (n,7) reaction 



rates in the literature sometimes differ from one another by 2-3 orders of magnitude. The 
calculated abundances of heavier elements based upon these rates can also differ from one 
another by 1-2 orders of magnitude. We carry out BBN calculations which include all of 
the recent compilations of reaction rates relevant to the production of isotopes including 
those that are heavier than 7 Li up to 18 (i.e. prito, Kajino, fc Uberhummer 1998|; |Mohr 



Hcrndl, & Oberhummcr 1999; Angulo et al. 19TJ9] ; |Herndl et al. 199"9" ; Wagcmans et ai 



1999| and any other previous estimates are considered). 
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3.5. Neutrino Degeneracy Parameters 

We have explored a broad range of the parameter space of neutrino-degenerate models. 
The main effects of the inclusion of either or v T degeneracy on BBN is an enhancement 
of energy density of the universe. The values for £„ and £„ T primarily affect the expansion 
rate. This means that £ v and £ Ut are roughly interchangeable as far as their effects on 
nucleosynthesis are concerned. Furthermore, we expect that the net total lepton number 
is small though the lepton number for individual species could be large. Hence, it is 
perhaps most plausible to assume that the absolute values of \£ v \ and |^ T | are nearly 
equal. Therefore, in what follows, we describe results for |£„ | = \£ Ut \ = . This reduces 
the parameter space to three quantities: Q^, £ Ve , and |£j,J = \£ Ut \- The calculations were 
conducted by first choosing a value for f2&/i| . It was then necessary to find a value of £„ e 
for which the light element constraints could be satisfied for some value of £^ T . 

3.6. Results 

As an illustration, Figure |5| shows a calculations of the helium abundance as a function 
of in a model with Qbhf = 0.3. It was found in this model that the helium constraint 
could be satisfied for £„ e in the range of 0.79 to 0.94. This figure is for £„ e = 0.9. Our 
results for £„ < 10 are consistent with those of Kang & Steigman (1992) if we run for a 
similar baryon-to-photon ratio r] w = 5 or 100 and utilize the neutron half life adopted in 
that paper. 

In our calculations for helium (Figure [5]) and other light nuclides (Figure ||) we see two 
different jumps in the nucleosynthesis yields. There is a small one for r w 10 and a large 
one for w 11.4. These two jumps are due to the fact that we consider two degenerate 
neutrino species which decouple at slightly different temperatures (due to the presence of 
muons to interact with v^). Therefore, one species (u T ) decouples above the QCD transition 
for a smaller value for . However, this jump has a smaller effect on nucleosynthesis 
because so many degrees of freedom are still carried by the other degenerate neutrino 
species (z^) which decouples later. 

The main difference between our results and Kang & Steigman (1992), however, is that 
the helium abundance is significantly changed at high values of £„ T . This results from our 
higher decoupling temperatures and therefore the lower relic neutrino temperatures during 
primordial nucleosynthesis (See Figure f|). This effect of varying T^/T 7 for large on BBN 
has not been adequately explored in the previous studies. 

Figure |6] illustrates the effects on the other light-element abundances for this particular 
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parameter set. This figure shows that for moderate values of £„ T the main effect is that 
weak reaction freeze-out occurs at a higher temperature. The resultant enhanced n/p ratio 
increases the abundances of the neutron-rich light elements, D, 3 H and 7 Li, while the 7 Be 
abundance decreases. 

Regarding 7 Li and 7 Be, the enhanced expansion rate from neutrino degeneracy affects 
the yield of A = 7 elements in two different ways. These elements are produced mainly 
by the nuclear electromagnetic capture reactions: t(a,7) 7 Li and 3 He(a,7) 7 Be. Hence, the 
production of these elements begins at a later time and lower temperature than the other 
light elements because they require time for a significant build up of the reacting A = 3 
and 4 elements. In a neutrino-degenerate universe, however, the increased expansion rate, 
hastens the freeze-out of these reactions from nuclear statistical equilibrium (NSE) resulting 
in reduced A = 7 yields relative to the nondegenerate case. However, the enhanced n/p 
ratio also increases the tritium abundance in NSE. This effect tends to offset the effect of 
rapid expansion on the production of 7 Li. The net result is more 7 Li production. 

As in the case of primordial helium, there is a rapid change in the nucleosynthesis 
yields on Figure |6] once the and v T decoupling temperatures separately exceed the epoch 
of the QCD phase transition. The ensuing lower neutrino temperature during primordial 
nucleosynthesis then resets the abundances to those of a lower effective degeneracy. The 
two discontinuities in Figure § at = 10 and 11.4 correspond to the points at which 
respectively the tau or muon neutrino decoupling temperatures separately exceed the QCD 
phase transition temperature (cf. Figure ^). 



Figure [7(a)| summarizes the allowed regions of the £ Ue vs. plane based upon the 
various indicated light element constraints in a universe with ^b^Jo = 0.1. The usually 



identified allowed region (cf. |Kohri, Kawasaki, fc Sato 1997[ ) for small £„ a ~ 0.2 and ~ 2 



is apparent. Indeed, it has been argued (cf. [Kohri, Kawasaki, fc Sato 1997] ) that such 
degeneracy may be essential to explain the differences in the constraints from primordial 
helium and deuterium. 



Figures |7( b )| and |7(c)| show the same plots, but for fl\,h\ Q = 0-2 and ^b^so = 0-3, 



respectively. The decline in the primordial deuterium abundance for models in which 
£ w > 10 allows for new regions of the parameter space in which the light element 
constraints can be accommodated. This new allowed region for large degeneracy persists as 
the baryon density is increased. 

Figure [8] highlights the basic result of this study. It shows that there exists a single 
parameter (either or J?b^Jo) family of neutrino degenerate models allowed by BBN. For 
low J^b^Jo models, only the usual low values for £„ e and £ u are allowed. Between O h hl 
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0.187 and 0.3, however, more than one allowed region emerges. For Q h hj > 0.4 only the 
large degeneracy solution is allowed. Neutrino degeneracy can even allow baryonic densities 
up to ^b^5o = 1- This result has been noted previously (cf. |Kang fc Steigman 1992| : 
Starkman 1992 ). What is different here is that the high i^b^lo models are made possible 
for smaller values of £„ e due to the higher neutrino decoupling temperatures deduced in 
the present work. As we shall see, this allows the circumvention of other cosmological 
constraints as well suggesting that baryons and degenerate neutrinos might provide a larger 
contribution to the universal closure density than has previously been derived. 

Figures PHTT| illustrate the elemental abundances obtained in the family of allowed 
models as a function of Ob^lo- Figure [D] in particular allows us to consider whether 
there exists an abundance signature in other elements which might distinguish this new 
degenerate neutrino solution from standard BBN. For the most part the yields of the light 
and heavy species are similar to those of the standard non-degenerate big bang. However, 
the boron abundance exhibits an increase with baryon density due to alpha captures on 7 Li. 
Thus, in principle, an anomalously high boron abundance together with beryllium and 6 Li 
similar to that expected from standard BBN might be a signature of neutrino-degeneracy. 



4. OTHER COSMOLOGICAL CONSTRAINTS 

We have seen that a new parameter space in the constraints from light elements on 
BBN emerges in neutrino-degenerate models just from the fact that the relic neutrino 
temperature is substantially diminished when degeneracy pushes neutrino decoupling to an 
earlier epoch. The viability of this solution requires large neutrino degeneracy. Hence, it 
becomes necessary to reexamine constraints posed from other cosmological considerations. 



4.1. Structure Formation 

It has been argued ( [Kang fc Steigman 1992] ) that large neutrino degeneracy is ruled 
out from the implied delay in galaxy formation in such a hot dark matter universe. This 
argument is summarized (Kang fc Steigman 1992) as follows: 



Neutrino degeneracy speeds up the expansion rate by a factor 

S 2 (6,) = 1 + 0.135(F(£„) - 3) « 0.135F(£,) , (12) 
where F(£ v ) is an effective energy density factor ( [Kang fc Steigman 1992| ) for neutrinos. 
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For massless neutrinos we have from Equations [3] and [5|, 

(13) 



L 7 



7 V 7T / 7 V 



7T 



where (T^/T^nuc is the neutrino-to-photon temperature ratio before e ± pair annihilation. 
This ratio is unity in the standard model with little or no degeneracy. 

With the increased mass-energy in degenerate neutrinos, the time of matter-radiation 
equality occurs at smaller redshift and there is less time for the subsequent growth of 
fluctuations. The redshift z eq (£ v ) for matter-radiation equality for a neutrino-degenerate 
universe can be written in terms of the redshift for a nondegenerate universe and the 
speed-up factor, 

1 + z eq &) = S " 2 (l + z eq {i v = 0)) , (14) 

For the present matter-dominated universe without neutrino degeneracy and Qm < 1 
we have 1 + z eq (£ u = 0) < 1.06 x 10 4 ( |Kang fc Steigman 1992| ). Furthermore, one 
demands that the fluctuation amplitude A(z) grows at least linearly with redshift, 
i.e. A(z) > (1 + z eq )A(z eq ). One also requires that the amplitude at least reaches unity by 
the present time, 

1 £ (1 + z eq ) < . (15) 



A(z eq ) A(z eqJ 
Thus, the requirement of sufficient growth in initial perturbations places a bound on the 
speed-up. Namely, 

1 06 x 10 4 

S 2 < — 7-t < 1.06 x 10 4 A(z eq ) . (16) 

J- ~r Z eq \<,v) 

Then requiring that A(z eq ) < 10~ 3 QStcigman fc Turner 1985|) leads to the constraint that 
F(f„) < 10.6/.135 w 79. 

Figure [1^ shows the F(£ u ) and F(£ Ui ) calculated as a function of £ u for the allowed 
models of Figure j|. Since our interesting parameter regions in Figure || satisfy £ Ue <^ ^ v , 
only v^ T contribute significantly to the total F(£ u ). Also shown are values of F(£ Vi ) if only 
one neutrino species was degenerate. The dashed line gives the constraint F(£ v ) < 79 from 
Kang & Steigman (1992). For comparison, the two-dot dashed line also shows the F(£ v ) for 
a single species from Kang & Steigman (1992). For low values of < 4, our F(£ Vi ) values 
for a single degenerate species are nearly identical to those of Kang & Steigman (1992). 
However, as £ Vi increases, our curves are lower due to the fact that we treat the annihilation 
epochs continuously (cf. Figure f|) rather than discretely, except for the QCD transition. 
By chance, our limit in the low degeneracy range for two degenerate neutrinos, (£ v T < 6.2) 
is close to the single species limit (^f" 5 ^ 6.9) of Kang & Steigman (1992). In the present 
work, however, the single species limit for moderate degeneracy increases to < 8.2. 
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Moreover, in the present work we now find that there also exists a new allowed region 
with 11.4 < T < 13 for which this growth constraint is satisfied even for two degenerate 
species. For only one degenerate species the new allowed range expands to 11.4 < £„ < 16. 
An upper limit of £„ r < 13 for two degenerate neutrino species corresponds to allowed 
BBN models with a baryon fraction as large as fif^io ^ 0.25 (cf. Figure ||). 



4.2. Neutrino Mass Constraint 

Figure O summarizes the neutrino contribution to the closure density Q u as a function 



of neutrino degeneracy and mass. This figure assumes the plausible model of nearly 
degenerate and v T masses and negligible v e mass. For this figure Q u refers to the 
combined contributions from both and v T . The contribution changes for large degeneracy 
due to the lower present-day neutrino temperature. This figure can be used to constrain 
the masses of the v and r neutrino types in different cosmological models. For example, if 
we assume a model with Qbh\ Q = 0.1, = 0.6, and Q u = 0.3, then we would find that 
the masses of both the and v T must be <2 eV, if these two species are to provide the 
neutrino contribution to the closure density. 



Studies of large scale structure also constrain neutrino masses and degeneracy in 
hot-dark-matter (HDM) and mixed-dark-matter (MDM) models. Indeed, at least some 
neutrino mass may presently be required to account for the observed power spectrum of 
galactic and microwave background structure. It has been argued ( |Primack et al. 1995| ) 



from considerations of structure formation in the early universe that two neutrino flavors 
(z^, z/ T ) may have a rest mass of 2.4 eV, compatible with all neutrino oscillation experiments. 
This postulate solves the main problem of cold dark matter (CDM) models, i.e. production 
of too much structure on small scales. Furthermore, Larsen & Madsen (1995) argue that 
neutrino degeneracy is required in a MDM model with 2.4 eV neutrinos to obtain an 
optimum fit to the power spectrum.. If we take 2.4 eV as the given mass of and z/ r , 
then for a Vt v < 0.9 ^lbh\ Q = 0.1 MDM cosmology we would deduce from Figure |l| that 
the maximum degeneracy for two species with this mass would correspond to £„ < 2.5 
similar to the values used in Larsen & Madsen (1995). We suggest that in light of the 
present results a similar study of the galactic power spectrum for £„ « 11 and m„ ~ 0.1 
eV should be undertaken as well. 
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4.3. Cosmic Microwave Background Constraint 

Perhaps, the most stringent remaining constraint on neutrino degeneracy comes from 
its effect on the cosmic microwave background. Several recent works ( [Kinney fc Kiotto 



1999|; |Lesgourgues fc Pastor 199~9| ; [Hannestad 2000|) have shown that neutrino degeneracy 



can dramatically alter the power spectrum of the CMB. The essence of this constraint is 
that degenerate neutrinos increase the energy density in radiation at the time of photon 
decoupling in addition to delaying the time of matter-radiation energy- density equality as 
discussed above. One effect of this is to increase the amplitude of the first acoustic peak in 
the CMB power spectrum at I « 200. For example, based upon a x 2 analysis ( [Lineweaverj 
|fc Barbosa 1998Q of 19 experimental points and window functions, Lesgourgues and Pastor 
(1999) concluded that £ v < 6 for a single degenerate neutrino species. 

However, in the existing CMB constraint calculations (Kinney and Riotto 1999: 
Lesgourgues and Pastor 1999; Hannestad 2000) only small degeneracy parameters with the 
standard relic neutrino temperatures were utilized in their derived constraint. Hence, the 
possible effects of a diminished relic neutrino temperature for large neutrino degeneracy 
need to be reconsidered. To investigate this we have done calculations of the CMB power 
spectrum, AT 2 = 1(1 + l)Ci/2n based upon the CMBFAST code of Seljak & Zaldarriago 
(1996). 

For the optimum neutrino-degenerate models ~ 11) and a neutrino contribution 

Q u < 0.25 we deduce from the solid curve on Figure ^| that the neutrino mass is m v < 0.1 
eV and therefore unimportant during the photon decoupling epoch. Therefore, we only 
consider massless neutrinos here. For massless neutrinos it can be proven ( [Lesgourgues fc| 
Pastor 1999) that the only effect of neutrino degeneracy on the CMB is to increase the 



background pressure and energy density of relativistic particles (cf. Eqs. 1-3). We have in 
this way explicitly modified CMBFAST code to account for the contribution of massless 
degenerate neutrinos with a relic temperature ratio y v = Tj/T u as given in Figure [|. 

We have evaluated \ 2 f° r fits to the CMB power spectrum, based upon the "radical 
compression" technique as described in Bond, Jaffe & Knox (2000). We have used the latest 
69 observational points and window functions available from the web page given in that 
paper. The advantage of this approach is that the non-Gaussian experimental uncertainties 
in the power spectrum are correctly weighted in the evaluation of the goodness of fit. 

For the purposes of the present study, we take as a benchmark the "All" case best 
fit 0=1 model of Dodelson & Knox (2000) who derived cosmological parameters based 
upon this same data set and compression technique. Although there is some degeneracy 
in the cosmological parameters they deduced an optimum fit to the power spectrum for 
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n b h 2 = 0.019, H = 65 km sec" 1 Mpc" 1 , fi A = 0.69, £l M = 0.31, r = 0.17, and n = 1.12, 
where Qm is the total matter contribution, r is the reionization parameter, and n is the 
"tilt" of the power spectrum. This benchmark is plotted as the dashed curve in Figure [TJ|. 
For this case we find \ 2 = 101. [Note that our \ 2 is slightly different from that quoted 
in Dodelson & Knox (2000) because we use different binning of the power spectrum]. For 
comparison the "radical compression" of the CMB data into 14 bins used in this work is 
also shown ( Pond et al. 2000] ). 



Rather than to do an exhaustive parameter search we have taken an approach similar 
to Lesgourgues & Pastor (1999). That is, we fix several representative cosmological 
models and then study their goodness of fit to the CMB data. The best case for large 
neutrino degeneracy will be for a value of the degeneracy parameter r such that neutrino 
decoupling occurs just before the QCD phase transition. This is the value for which the 



relic neutrino energy density is a local minimum (cf. Figure [13]). For the present work this 



corresponds to £„ = 11.4, £„ e = 0.73, fib^io — 0.187 models. In what follows we fix £ u , 
£„ e , and ^ b /i| at these values and refer to this as the large degeneracy model. 

We have found [as did Lesgourgues & Pastor (1999)] that the currently favored 
= 0.7 models give a poor fit to the data even with no degeneracy. Adding neutrino 
degeneracy to an Q\ = 0.7 model only makes the fit worse. The main problem is that the 
first acoustic peak increases in amplitude and moves to larger /. Hence, even though a local 
minimum develops for large degeneracy, the \ 2 is substantially increased and large neutrino 
degeneracy is probably ruled out for = 0.7 models. 

For smaller VL\ a local minimum develops in the x 2 f° r both small values of degeneracy 
£„ « 1 and large degeneracy = 11.4. As pointed out in Lesgourgues & Pastor (1999), 
the best case for neutrino degeneracy is with Q\ = models. However, those models are 
probably ruled out by observations of type la supernovae at high redshift ( |Garnavich, et 



|al. 1998| ; jPerlmutter et al. 1998a| ; Perlmutter et al. 1998b] ; |Riess et al. T99p . At the 



3cr confidence level for Q = 1 models, the type la Supernova data are consistent with 
^a = 0.7 ± 0.3 Hence, we take ^a = 0.4 as a plausible cosmological model which is 
marginally consistent with the type la results. Nevertheless, for purposes of illustration, 
we have also made a search for optimum parameters for matter dominated Q\ = 0, Q — 1 
models. 

The reason low Q\ models are favored is that they shift the first acoustic peak back 
to lower I. Larger values of H also slightly improve the fit by shifting the first acoustic 
peak to lower / and decreasing the baryon density for fixed ^ b /?| which lowers the peak 
amplitude. We take H Q = 65 ± 10 (h 50 = 1.3 ± 0.2) as a reasonable range ( Podelson fc| 
Knox 20"00| ), and therefore utilize h§o = 1.5 as the optimum Hubble parameter for the 
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neutrino-degenerate models. This implies that = 0.084 for the large degeneracy models. 
The ionization parameter does not particularly help the fits as it mainly serves to decrease 
the amplitude of both the first and second peaks in the power spectrum. We therefore set 
t = for the large degeneracy models. The only remaining adjustable parameter of the fits 
is then the tilt parameter n. Values of n slightly below unity also help with the amplitude 
and location of the first acoustic peak. 

The solid line on Figure [14] shows a = 0.4 model for which n = 0.78. For this fit 
Ax 2 = 27 which makes this large degeneracy model marginally consistent with the data 



at a level of 5.2a. The dotted line in Figure 111 shows the matter dominated ^a = best 



fit model with n = 0.83. For this fit A% 2 = 9 which makes this large degeneracy model 
consistent with the data at the level of 3a. 



As can be seen from Figure [14], a model with large neutrino degeneracy seems 
marginally acceptable based upon the presently uncertain power spectrum. The main 
differences in the fits between the large degeneracy models and our adopted benchmark 
model are that the first peak is shifted to slightly higher I values and the second peak is 
somewhat suppressed. It thus becomes important to quantify the amplitude of the second 
peak in order to constrain the large degeneracy models proposed herein. 

Indeed, after the present fits were completed a suppression of the second acoustic peak 
in the power spectrum was reported in the high-resolution BOOMERANG (Pernardls, P. et| 
al. 20001 ; |Lange et al."2000|) and MAXIMA- 1 QHanany et al. 200q : |Balbi et al. 20001 ) results. 



We have not yet analyzed the goodness of fit to these data as the experimental window 
functions are not yet available. In a subsequent paper we will examine the implications of 
those data in detail. 



For purposes of illustration, however, we compare the fit models of Figure [14| with 



the published BOOMERANG and MAXIMA- 1 power spectra in Figure |T|. Here one can 
clearly see that the suppression of the second acoustic peak in the observed power spectrum 
is consistent with our derived neutrino-degenerate models. In particular, the MAXIMA- 1 
results are in very good agreement with the predictions of the neutrino-degenerate 
cosmological models described herein. There is, however, a calibration uncertainty between 



these two sets ([Hanany et al. 200U|). If one only considers the BOOMERANG results alone, 



the diminished amplitude of the first acoustic peak probably tightens the constraint for 
low neutrino degeneracy models (cf. Hannestad 2000|) although even for this set alone, 
a high degeneracy model is probably still acceptable ( [Lesgourgues fc Peloso 20001) . It 
is clear, that these new data sets will substantially improve the goodness of fit for the 
neutrino-degenerate models ( [Lesgourgues fc Peloso 20001) . Moreover, both data sets seem 



to require an increase in the baryonic contribution to the closure density as allowed in our 
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neutrino-degenerate models. 

5. CONCLUSIONS 

We have discussed how the relic neutrino temperature is substantially diminished 
in cosmological models with a large neutrino degeneracy. We have shown that all of 
the BBN light-element abundance constraints (assuming some destruction of 7 Li) can 
be satisfied for a single-parameter family of cosmological models in which significant 
neutrino degeneracies and large values of Vt^h 2 can exist. The requirement that large 
scale structure become nonlinear in sufficient time can also be satisfied for models with 
either moderate degeneracy < 6.2 and J^b^o ~ 0.22) or large neutrino degeneracy 

(11.4 < £ v < 13 and 0.187 < £2hhl < 0.25). We have also shown that even the constraint 
from neutrino-degeneracy effects on fluctuations of the cosmic microwave background 
temperature may be marginally avoided for models with < 0.4, £ v 11, and 
Q h hl Q « 0.2. 

At present, the power spectrum of the CMB is the most stringent constraint. 
Nevertheless, neutrino-degenerate models can be found which are marginally consistent at 
the 3-5a level. This tight constraint is due, at least in part, to a suppression of the second 
acoustic peak in the spectrum. It is therefore encouraging that the recent BOOMERANG 
and MAXIMA- 1 results suggest that such a suppression in the second acoustic peak may 
indeed occur in agreement with the expectations of the large neutrino degeneracy, high fl h 
models proposed here. 

Thus, high resolution microwave background observations become even more important 
as a means to quantify the limits to (or existence of) possible cosmological neutrino 
degeneracy. Based upon the current analysis, we conclude that all of the cosmological 
constraints on large neutrino degeneracy are marginally satisfied when a careful accounting 
of the neutrino decoupling and relic neutrino temperature is made. It will, therefore, be 
most interesting to see what further constraints can be placed on this possibility from the 
soon to be available space-based high resolution CMB observations such as the NASA MAP 
and ESA Planck missions. 
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Fig. 1. — Decoupling temperature To (in MeV) for v e (solid line), (dot-dashed line), 
and v T (dashed line) as a function of degeneracy parameter £ vv compared with the previous 
estimate (two-dot dashed line) of Kang & Steigman (1992) 



-28- 




Fig. 2. — Ratio of the energy density of muons (or any massive particle //) to photons as a 
function of the ratio of temperature to rest mass. This shows that even at a temperature of 
only ~ 20% of the rest mass, a significant fraction (~ 10%) of the energy density still resides 
in particle/antiparticle pairs. 
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Fig. 3. — An illustration of the effective degrees of freedom g e ff as a function of temperature 
for a standard big bang with 3 nondegenerate neutrinos. For this illustration we have 
assumed that the neutrinos do not decouple until T w MeV. The discontinuity at T = 150 
MeV is due to the QCD phase transition. 
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Fig. 4. — Ratio of relic neutrino temperature to photon temperature as a function of 
degeneracy parameter for each neutrino species for v e (solid line), (dot-dashed line), 
v T (dashed line), compared with the previous estimate (two-dot dashed line) of Kang & 
Steigman (1992) 




Fig. 5. — Helium mass fraction Y p as a function of degeneracy parameter £„ MT for O b /i§ = 0.3 
and £„ e = 0.9. Horizontal lines show the adopted observational constraints. 




Fig. 6. — The predicted abundances of the light elements as a function of the neutrino 
degeneracy ^ for f2 h hl = 0.3 and £ Ve = 0.9. 
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Fig. 7(a). — Contours of allowed values in the £„ e — plane for f2 b /i§ = 0.1, based upon 
the various light-element abundance constraints as indicated. The hatched region depicts 
the allowed parameters consistent with all light element constraints for this value of fl\yh\ Q . 
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Fig. 8. — Allowed values of £„ e and £„ for which the constraints from light element 
abundances are satisfied for values of ^bhl = 0.075, 0.1, 0.2, 0.3 and 1.0 as indicated. 
For large values of ^b^so > 0-3 the only allowed regions are for the large values > 20. 
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Fig. 9. — The predicted Helium abundance for allowed neutrino-degenerate models as a 
function of flhh\ - The v& l ues of £ Ve and £ u were taken from the central value in the 
allowed region determined by Obh^ i n Fig- H- Note that Y p is sensitive to the choice of £ Ve in 
the allowed region. 
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Fig. 10. — The predicted D/H, 3 He/H, and total A = 7 and 7 Li abundances for allowed 
neutrino-degenerate models as a function of J^b^Jo- The values of £ Ue and £„ were taken 
from the central value in the allowed region determined by ^b^so m Fig- H- 
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Fig. 11. — The predicted 6 Li, 9 Be, n B and A > 12 heavier element abundances as a function 
of J?b/?| . The values of £„ e and ^ were taken from the central value in the allowed region 
determined by ^b^lo m Fig- ||- 



-40- 




Fig. 12. — Contours of equal present energy density of massive degenerate neutrinos as a 
function of the degeneracy £ u and neutrino mass m u for Vt v = 0.1, 0.3, 0.5, and 1.0 as 
indicated. Each curve corresponds to different value of fl u as indicated. Shaded regions 
depict the allowed range of degeneracy for the two indicated values of f2b^lo = 0.1, and 0.2 
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Fig. 13. — Calculated neutrino energy density factors F(£ Ui ) as a function of degeneracy 
parameter for the three neutrino species. The dotted and dot-dashed curves display 
respectively F(£ v ) and F(£ Ut ) for the cases in which only one neutrino species is degenerate. 
Since our interesting parameter regions in Figure |] satisfy ^ Ve <C £ v , only v^ T contribute 
significantly to the total F(£ u ) (solid curve). The dashed horizontal line indicates the 
constraint on neutrino degeneracy from the requirement that sufficient structure develops 
by the present time. We also show the previous estimate (two-dot dashed curve) of Kang & 
Steigman (1992). 
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Fig. 14. — CMB power spectrum compared with calculated Q — 1 models. The points show 
the binning of 69 experimental measurements based upon the radical compression method of 
Bond et al. (2000). The dashed line shows the optimum model of Dodelson & Knox (2000). 
The solid line shows an Q\ = 0.4 model with three degenerate neutrinos £ v = 11.4, 
£j, e = 0.73 and flhh\ Q = 0.187 as described in the text. The dotted line is for an Q\ = 
model with the same degeneracy and baryonic parameters. 
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Fig. 15.— CMB power spectrum from the recent MAXIMA-1 (circles) and BOOMERANG 
(squares) binned data compared with calculated Q = 1 models of Figure 14. Note that the 
suppression of the second acoustic peak in the data is consistent with that predicted by the 
neutrino-degenerate models. 



